Transition from damage to fragmentation in collision of solids 
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We investigate fracture and fragmentation of solids due to impact at low energies using a two- 
dimensional dynamical model of granular solids. Simulating collisions of two solid discs we show 
that, depending on the initial energy, the outcome of a collision process can be classified into two 
states: a damaged and a fragmented state with a sharp transition in between. We give numerical 
evidence that the transition point between the two states behaves as a critical point, and we discuss 
the possible mechanism of the transition. 
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I. INTRODUCTION 



Fragmentation, i.e. the breaking of particulate ma- 
terials into smaller pieces is a ubiquitous process that 
underlies many natural phenomena and industrial pro- 
cesses. The length scales involved in it range from the 
collisional evolution of asteroids through the scale of ge- 
ological phenomena down to the breakup of heavy nu- 
clei [0J^]. In most of the realizations of fragmentation 
processes the energy is imparted to the system by im- 
pact, i.e. typical experimental situations are shooting 
a projectile into a solid block, free fall impact with a 
massive plate and collision of particles of the same size 
P-tnl. The most striking observation about fragmenta- 
tion is that the size distribution of fragments shows power 
law behavior independent on the microscopic interactions 
and on the relevant length scales, e.g. the charge distri- 
bution of small nuclei resulted from collisions of heavy 
ions exhibits the same power law behavior as the size 
distribution of asteroids Experiments revealed 

that the power law behavior of fragment sizes is valid 
for a broad interval of the imparted energy |^[ll|l , which 
was also reproduced by computer simulation of sophis- 
ticated microscopic models ]T^-|T7||. The observation of 
power law size distributions without control parameter 
initiated the idea of self-organized criticality p|- pd] , p^| , ^9[ 
in fragmentation and gave rise to numerous studies to 
understand the dynamic origin of the power law [p2|-p5[ . 
Hence, during the past years much information has been 
accumulated about fragmentation processes in the im- 
parted energy range where a power law size distribution 
occurs but the limiting case of low energies is still not 
explored. 



Beside the general interest in fragmentation processes 
one can also mention other fields where fracture and frag- 
mentation of solid particles due to impact play an impor- 
tant role. It is well known that in the flow of granular 
materials a large part of the kinetic energy of the grains is 
dissipated in the vicinity of their contact zone during the 
collisions. Beside the viscous and plastic effects, the dissi- 
pation by damaging is also an important source of energy 
loss in the flow. Collision of particles occurs also in the 
solar system in planetary rings. In this case the energy 
dissipation due to impact damage might also influence 
the large scale structure formation in the rings |2q| . On 
larger length scale in the solar system, the so-called colli- 
sional evolution of asteroids due to subsequent collisions, 
and the formation of rubble piles in the asteroid belt is 
still a challenging problem || . Among industrial applica- 
tions the breakup of agglomerates in chemical processes 
can be mentioned. Due to experimental difficulties, the 
computer simulation of microscopic models is an indis- 
pensable tool in the study of these impact phenomena 

6M1- 

In the present paper we want to elaborate the impact 
fracture and fragmentation of solids at low imparted en- 
ergy using a two-dimensional dynamical model of break- 
able granular solids. Simulating collisions of two solid 
discs we show that, depending on the imparted energy, 
the outcome of a collision process can be classified into 
two states: a damaged and a fragmented state with a 
sharp transition in between. Analyzing the energetics of 
the impact and the resulted fragment size distributions 
we give numerical evidence that the transition point be- 
tween the damaged and fragmented states behaves as 
a critical point. The transition proved to be the lower 
bound for the occurrence of power law size distributions. 
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The possible mechanism of the transition between the 
two states is discussed. In spite of the specific features of 
the system studied here, most of our results can be con- 
sidered generally valid for impact phenomena mentioned 
above. 

After giving a short summary of the main ingredients 
of our model in Sec. II., the numerical results concerning 
to the energetics of the collision process and to the size 
distribution of fragments will be presented in Sees. III. 
and IV., respectively. In Sec. V. we discuss the possi- 
ble mechanism of the transition and some general con- 
sequences of our work for other types of fragmentation 
phenomena. 



II. THE MODEL 



gons until the entire system relaxes, i.e. there is no break- 
ing of the beams during some hundreds consecutive time 
steps and there is no energy stored in deformation. For 
more details of the model's definition see Refs. ]T^ , p^ |. 

We have applied the model to study shock fragmenta- 
tion of solids in various experimental situations. Namely, 
simulations were performed to study the fragmentation of 
a solid disc caused by an explosion in the middle |l2| , the 
breaking of a rectangular block due to the impact with a 
projectile [12]], and the collision of two macroscopic bod- 
ies (discs) O. The model proved to be successful in re- 
producing the experimentally observed subtleties of frag- 
menting systems, e.g. the power law mass distribution 
of fragments was found independent on the initial con- 
ditions, with an exponent in the vicinity of two, slightly 
depending on the initial energy |12jjR 



Recently, we have worked out a two-dimensional dy- 
namical model of deformable, breakable granular solids, 
which enables us to perform molecular dynamics simu- 
lation of fracture and fragmentation of solids in various 
experimental situations [12 13 1. Our model is an exten- 



sion of those models which are used to study the behavior 
of granular materials applying randomly shaped convex 
polygons to describe grains [ p7fl . To capture the elas- 
tic behavior of solids we connect the unbreakable, unde- 
formable polygons (grains) by elastic beams. The beams, 
modeling cohesive forces between grains, can be broken 
according to a physical breaking rule, which takes into 
account the stretching and bending of the connections. 
The breaking rule contains two parameters t e , te control- 
ling the relative importance of the stretching and bend- 
ing breaking modes, respectively. The energy stored in a 
beam just before breaking is released in the breakage giv- 
ing rise to energy dissipation. The average value of the 
energy dissipated by the breakup of one contact defines 
the crack surface energy E s in the model solid as a func- 
tion of the breaking parameters t e , and tg>. At the broken 
beams along the surface of the polygons cracks are gen- 
erated inside the solid and as a result of the successive 
beam breaking the solid falls apart. The fragments are 
defined as sets of polygons connected by remaining intact 
beams. In the framework of this discrete model one can 
introduce the notion of the binding energy Eh as the en- 
ergy released in the case of complete disintegration, i.e. 
in the case when all the fragments are single polygons. 
Hence, 



E h = Nr. ■ E K 



(1) 



where N c denotes the total number of grain-grain con- 
tacts in the solid. Note that Eh is proportional to the 
volume of the sample and can be considered as a natural 
energy scale of the system. 

The time evolution of the fragmenting solid is obtained 
by solving the equations of motion of the individual poly- 



III. DAMAGE AND FRAGMENTATION 

In the present paper we apply our model to explore 
the properties of impact fragmentation processes at low 
imparted energy. For this purpose we carried out molecu- 
lar dynamics simulation of the collision of two solid discs 
of the same size, similar to Ref. (juj). The disc-shaped 
granular solid was obtained starting from the Voronoi- 
tessellation of a square [§8|,^9) and cutting out a circular 
window in the middle. This way of construction gives 
rise to a certain surface roughness of the particles. The 
schematic representation of the collision of two particles 
is depicted in Fig. [l| For simplicity, in the present studies 
only central collisions were considered, i.e. the value of 
the impact parameter b was set to zero in all the simula- 
tions, and only the size of the particles R and the surface 
energy E s were varied. In the simulations the energy- 
release by beam breaking is the only possible source of 
dissipation. The energy of the collision Eq is defined as 
the total initial kinetic energy of the colliding bodies 



Eq = R 2 npvQ, 



(2) 



where p is the mass density and vq denotes the initial 
velocity of the particles. When Eq is smaller than the 
so-called damage threshold Ed, the collision results sim- 
ply in elastic rebound without internal damage (without 
breaking of beams). To achieve damaging, Eq has to sur- 
pass the damage threshold Ed, which is determined by 
the surface energy E s of the model. 

It is generally accepted that impact fragmentation phe- 
nomena exhibit so-called energy scaling, i.e. the result 
of the fragmentation process only depends on the value 
of the specific energy defined as the imparted energy di- 
vided by the total mass of the system. To characterize 
the collision events we introduce a dimensionless param- 
eter r] with the definition 
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This choice of parameter r\ has the advantage with re- 
spect to the specific energy that besides the dependence 
on the global size it includes also the specific material 
type through the surface energy E s . Since in the present 
study we focus on the evolution of damage with increas- 
ing impact energy, simulations were performed at fixed 
values of R and E s varying 77 in the range 770 < 77 < 1.5. 

Here 770 corresponds to the damage threshold 770 — \J~%^- 
Our results concerning to impact fragmentation at 77 > 
1.5 can be found in Ref. jli| including detailed analysis 
of the dynamics of the collision process and comparison 
to experiments. The values of the most important pa- 
rameters of the present simulations are summarized in 
Table |. 

In Fig. H we present the final breaking scenarios (re- 
laxed states) of the collision process obtained by simula- 
tions at four different values of 77. It can be observed that 
for 77 ~ r/o (Fig. ||a) cracking by beam breaking mainly 
occurs in the vicinity of the contact surface of the two 
bodies and the bulk remains practically intact. Since 
this gentle collision does not cause size reduction of the 
two particles, this case can be considered as an inelastic 
impact of discs, where energy dissipation is solely due 
to internal damage. Increasing the impact energy Eq, 
i.e. increasing 77, gives rise to more broken beams in Fig. 

and the solids break into pieces. Around the impact 
site of the bodies the fragments are smaller (single poly- 
gons, some pairs and triplets), and there are a few much 
larger fragments the size of which is still comparable to 
the original size of the bodies. Further increase of 77 (Fig. 

, d) mainly results in breakup of the large fragments 
into smaller ones giving rise to fragmentation of the en- 
tire solids. Based on the above qualitative picture, we 
classify the outcome of a collision into two states, i.e. 
damaged and fragmented states distinguished by the size 
of the largest remaining piece. The colliding bodies are 
considered to be damaged when the size of the largest 
piece is comparable to the original size of the bodies (see 
Fig. ||a, 6), while the fragmented state is characterized by 
the absence of such large pieces (Fig. ||c, 0?). In the fol- 
lowing we point out that the above qualitative picture is 
also supported by the quantitative analyzis of the ener- 
getics of the collision process and of the behaviour of the 
resulting size distribution of fragments. Namely, evalu- 
ating the energy released by the breaking of beams, and 
quantities related to the size distribution of fragments we 
will justify that the behaviour of the colliding system has 
two different regimes with a sharp transition in between. 

The energetics of the collision process corresponding 
to the samples of Fig. || is summarized on Figs. ||, and 
|j. Fig. U shows the energy En released by beam break- 
ing, together with the remaining kinetic energy Eq — En 
stored in the motion of fragments as a function of 77. For 



the purpose of comparison, the surface energy E s of the 
model, the binding energy of the sample E^ and the en- 
ergy of the impact E are also plotted in the figure. When 
the impact energy Eq surpasses the damage threshold Ed 
(77 surpasses 770) the released energy En takes its small- 
est value, which is equal to the surface energy E s of the 
model (only one beam is broken). The largest possible 
value of En is equal to the binding energy of the sample 
Eb, which can only be reached in the case of complete 
disintegration. According to the simulations, this limit- 
ing case is hard to achieve, since on the surface of the 
colliding discs opposite to the impact site, small frag- 
ments comprising a few polygons can always escape. It 
can also be seen in the figure that the impact energy Eq 
producing complete disintegration has to be much larger 
than the binding energy E\,. Note that similar qualita- 
tive behaviour was found simulating free fall impact of a 
disc with a massive plate |l6|] . 

The most interesting observation of Fig. ||| is that the 
curve of En versus 77 is composed of two distinct parts, 
i.e. a rapidly increasing low energy part, and a slowly 
increasing high energy part. These two regimes are sep- 
arated by a point where En is practically equal to the 
remaining kinetic energy Eq — En, which implies that at 
this point half of the total kinetic energy Eo is released 
by cracking. The two regimes can be better observed in 
Fig. |] where the ratio e n of En and Eq is plotted against 
77. The curve of en has a maximum, the value of which is 
in the vicinity of 0.5. This maximum separates a rapidly 
increasing and a slowly decreasing part, which we asso- 
ciate to the damaged and fragmented regimes introduced 
above, respectively. We identify the transition point be- 
tween the two states with the position of the maximum 
called fragmentation threshold r\ c . 

The separation of the two regimes and the identifi- 
cation of the transition point is also supported by the 
behaviour of the mass of the largest fragment as a func- 
tion of 77. In Fig. P we present the sum of the mass of 
the largest fragments of the two bodies M max — M^ ax + 
M max normalized by the total mass M tot = M A + M B 
versus 7/. The monotonically decreasing function has a 
distinct curvature change, the position of which r\ c — 0.36 
coincides with the maximum of the energy release curve 
en. This also confirms that the behaviour of the system 
sharply changes at a specific value of 77 that we called 
fragmentation threshold ?7 C . 

To clarify the influence of the overall size of the collid- 
ing bodies on the behaviour of the fragmenting system, 
simulations were performed varying the radius R of the 
particles between 10 and 30 cm. Note that the average 
size of the randomly shaped unbreakable polygons is 1cm. 
Fig. ^ shows en a s a function of 77 for several values of 
R. It can be observed that the qualitative behaviour of 
the curves is universal, it is independent on the system 
size, but the position of the transition point r\ c is slighly 
shifted downward with increasing R. 
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IV. THE SIZE DISTRIBUTION OF FRAGMENTS 

To reveal the nature of the transition from the dam- 
aged to the fragmented state, the evolution of the frag- 
ment size distribution with varying impact energy is a 
crucial point. Since the dissipated energy Er is propor- 
tional to the total surface of cracks it is expected that 
the two regimes should clearly show up also in the type 
of disintegration of the solids, i.e. in the behaviour of the 
size distribution of fragments. In the present chapter we 
demonstrate that in the vicinity of rj c the behaviour of 
the fragment size distribution shows strong similarities 
to that of the size distribution of clusters in systems un- 
dergoing second order phase transition. As an example, 
one can mention the behaviour of the cluster size distri- 
bution in percolation around the critical point p c , or the 
size distribution of droplets in the vicinity of the critical 
temperature T c in the case of the liquid-gas phase tran- 
sition. Throughout this chapter we refer to percolation 
on finite lattices for the purpose of comparison. 

The fragment mass histograms F(m; rj) corresponding 
to the colliding system of the preceding chapter are shown 
in Figs. 0a, b for several values of rj below and above the 
transition point rj c . Here F(m; rf) denotes the number of 
fragments with mass m divided by the total number of 
fragments, averaged over 10 collision events having the 
same rj. In order to obtain the same statistics of the dis- 
tribution at all sizes a logarithmic bining was used, i.e. 
the bining is equidistant on logarithmic scale. The his- 
tograms have two cutoffs. The lower one is due to the 
existence of single unbreakable polygons while the upper 
one is determined by the finite size of the bodies. 

In Fig. H we have shown that a collision at rj just above 
rjo results in only a few small fragments and two big ones 
(with almost the same size) but no fragments are gen- 
erated in the intermediate mass region. Hence, in Fig. 
0a the corresponding mass distribution F(m; rf) has two 
peaks at small and at large mass with a large gap in be- 
tween. Increasing rj in Fig. 06, the main difference com- 
pared to Fig. 0a is that the peak of the large fragments 
gradually disappears, the large pieces break into smaller 
ones giving rise to fragments in the intermediate mass 
region. For rj — > rj c the mass distribution F(m; rf) in the 
intermediate region tends to a power law: 



F(m) ~ m T . 



(4) 



The value of the exponent obtained at rj = rj c is r = 



2.27±0.05. Our former simulations showed (see Ref. |13|) 
that the power law behavior of F(m; rf) remains valid for 
•q > rjc, but the value of the exponent r is a decreasing 
function of 77, in agreement with most of the experimental 
observations pH7|. 

To clarify the effect of the overall size of the colliding 
bodies R on the shape of the mass distribution F(m; rf) 
and on the value of the exponent r, in Fig. @ we plotted 



the fragment size distributions at the transition point rj c 
for three different values of R. It can be seen that in- 
creasing R the power law region of F(m; rj c ) gets wider 
and the relative height of the hump at large fragments de- 
creases, but the value of the exponent r does not change. 
It is important to note that in the case of percolation on 
finite lattices the size distribution of clusters shows the 
same shape and dependence on the system size at the 
critical point fttj]. 

More insight into the evolution of the shape of F(m; rj) 
can be obtained by studying the moments of the distri- 
bution as a function of rj. The fcth moment Mk(rf) of the 
histogram F(m; rf) is defined as 



Mk{rj) =Y^m k ■ F(m;rf). 



(5) 



In the case of critical phenomena like percolation or the 
liquid-gas transition the moments Mk of the cluster size 
distribution with k > 1 diverge at the critical point in 
the thermodynamic limit 



Mi. 



->ik 



(6) 



where e denotes the distance from the critical point, i.e. 
e = p — Pc for percolation or e = T — T c for the liquid-gas 
transition. In a finite system the moments Mk have a 
finite maximum at the transition point. Assuming gap 
scaling for F(m) |3(J 



F(m) ~ m T f(m a e) 



(7) 



the moment exponents jik can be expressed in terms of 
r and a 



1 



(8) 



and the behavior of the system in the vicinity of the tran- 
sition point can be characterized by only two independent 
exponents t and a. 

To test whether our small systems exhibit some trace 
of this behavior we evaluated so-called single event mo- 



ments Ml defined as 



Ml = 'm k 



nP (m) 



(9) 



where the upper script j refers to the jth collision event, 
n J (m) denotes the number of fragments with mass m in 
event j and the prime indicates that the sum runs over 
all the fragments excluding the largest ones of each of the 
two colliding bodies. The study of single event moments 
was first suggested by X. Campi to reveal properties of 
fragmentation processes with varying imparted energy 
pl|-p3|. As it was demonstrated in Refs. pl|-p3[, the 
usage of single event moments defined by Eq. (P) instead 
of Eq. (^) has the advantage that part of the analysis, 
e.g. the study of the correlation of moments in scatter 
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plots, can be carried out without ordering the collision 
events according to a parameter like in our case rj. Fol- 
lowing Campi's ideas, we evaluated the ratio of M% and 



M{ 



Ml 

M{ £Lm-»*(m) 



(10) 



which is equal to the average fragment size M in event 
j. In Fig. | (M 3 2 /M{\ is plotted as a function of r\ for 
several different values of the system size R. The brack- 
ets < . . . > denote that each data point was obtained 
as an average over 10 events having the same r/. One 
can observe that (mI/M(\ has a distinct maximum the 
position of which coincides with the maximum of the en- 
ergy release curve in Fig. |^, within the precision of 
the calculations. It is very important to note that in- 
creasing the system size R the peak of (^M^/M(^ gets 

sharper, i.e. the height of the peak increases while the 
width of the peak decreases, which is also typical for a 
critical point like the percolation threshold occurring in 
finite size systems pO| ) Pl| , p^| . 

The validity of gap scaling Eq. (Q) for the mass dis- 
tribution F(m; rf) can be easily tested by making scatter 
plots of moments M k , i. e. by plotting M k against M° k , 
with k, k > 1, then checking the correlation of the mo- 
ments and the validity of the relation of the exponents 
Eq. (^|). The scatter plot of two pairs of moments is 
For the purpose of normalization the 
> 1 are divided by M{ . In Fig. [lO] the 



shown in Fig. 10 
moments with k 



collision events close to the transition point r\ c are repre- 
sented by points with the largest values of the moments. 
Events with r\ far below or far above r/ c result in smaller 
values of M k , hence, the corresponding points in Fig. 
|Io| fall closer to the origin. One of the most remarkable 
features of Fig. |l^ is that the moments are strongly cor- 
related in the region corresponding to the vicinity of the 
transition point rj c and a power law seems to be a reason- 
able fit to the scattered points. By this fitting procedure 
one obtains numerically the ratio of the moment expo- 
nents (ik and (i k < since 



k 



k,k > 1. 



(11) 



The value of the exponents obtained as the slope of 
the straight lines in Fig. [l(] are ^3/^2 — 2.28 ± 0.1 
and fJ-s/n2 — 4.95 ± 0.17, which are in agreement with 
the corresponding predictions 2.36 and 5.11 substituting 
t = 2.27 into Eq. (g). 

Hence, we conclude that, within the limited precision 
available in our small system, the mass distribution of 
fragments fulfills the gap scaling relation Eq. (Q), and the 
damage-fragmentation transition occurs as a continuous 
phase transition. The control parameter of the transition 



can be identified with the parameter 77, and the quantity 
characterizing the size reduction M max / M to t can be con- 
sidered as the order parameter. To complete the char- 
acterization of our system around r/ c we determined nu- 
merically the value of the order parameter exponent (3 
defined as |35|l 



M„ 



M tl 



\V-Vcf, 



V < Vc 



(12) 




and the value of the exponent 7 characterizing the aver- 
age fragment size 



(13) 



The power law behavior given by Eq. (|l|) is valid 011 
both sides of the critical point r\ c . (Note that 7 = /i2-) 
Using the data presented in Fig. p|we plotted M max /M to t 
versus \rj — r] c \ for r\ < r\ c in Fig. FLU It can be seen that a 
power law with an exponent j3 — 0.11 ± 0.02 gives a rea- 
sonable fit to the data. In Fig. (^M|/M^ belonging 
to the system size R = 20cm of Fig. ^ is depicted as a 
function of | rj — r\ c | for ?y > rj c . Again a power law with an 
exponent 7 = 0.26 ± 0.02 can be fitted with a reasonable 
quality. Since only two of the critical exponents are inde- 
pendent, we can check the consistency of our description 
by pairing the exponents r, /3, and 7 and using the well 
known scaling laws [ j30|]35|l of critical phenomena to ob- 
tain the value of the other exponents. It can be observed 
in Table [n] that the numerical values of the exponents 
obtained by starting from the three possible pairs agree 
well with each other within the error bars, which gives 
further justification of our phase transition picture. 



V. DISCUSSION 

In the present paper we applied our dynamical model of 
granular solids to study impact fragmentation phenom- 
ena at low values of the imparted energy. Analyzing the 
energetics of the fragmentation process and the resulting 
size distribution of fragments we identified two distinct 
final states of the impact process, i.e. damaged and frag- 
mented states with a sharp transition in between. With 
a detailed study of the behavior of the fragment mass 
distribution in the vicinity of the transition point and its 
dependence on the finite particle size we gave numerical 
evidence that the transition point behaves as a critical 
point and the damage-fragmentation transition occurs as 
a continuous phase transition. The control parameter of 
the transition was chosen to be the dimensionless ratio 
77 of the energy of impact and the binding energy of the 
sample, and the order parameter was associated to the 
mass of the largest fragment divided by the total mass. 
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Based on the time evolution of the collision process ob- 
tained by the simulations, the following qualitative pic- 
ture of the mechanism of the transition can be estab- 
lished: when two macroscopic solids collide damage first 
occurs in the surroundings of the contact zone of the 
bodies. After this contact damage a compressive elas- 
tic pulse expands radially through the bodies the ampli- 
tude of which strongly depends on the amount of primary 
damage occurred in the contact zone. The pulse reflects 
back from the free boundary and after reflection tensile 
forces arise in the solids resulting in cracks in the bulk. 
The amount of breaking in the bulk is mainly determined 
by the amplitude of the initial pulse remaining after con- 
tact damage. The fragmentation of the colliding solids 
occurs above a specific value of the impact energy where 
the amplitude of the remaining pulse will be sufficient to 
give rise to the complete break up of the bulk into pieces. 
At this specific energy value the behavior of the colliding 
system sharply changes as we demonstrated by analyzing 
the energetics and the size distribution of fragments. 

The idea of the existence of a so-called fragmentation 
phase transition was first addressed in the field of nu- 
clear physics where the fragmentation of heavy nuclei 
due to impact is extensively studied (for a recent review 
see Ref. |36|]). Using percolation based ideas it was shown 
that the disassembly of excited nuclei possesses a contin- 
uous phase transition when the imparted energy is varied 
pl| p4| . The order parameter of the transition was asso- 
ciated to the size (charge) of the largest fragment. 

Recently, the penetration of a steel ball into a solid 
plate with varying impact energy was studied experi- 



mentally |37|. It was reported that the so-called ballistic 
limit, i.e. the impact energy where the perforation of the 
plate occurs, behaves as a critical point, and the perfo- 
ration occurs as a continuous phase transition the order 
parameter of which was chosen to be the mass of ejecta 
expelled from behind the target. It was argued that the 
mechanism of the phase transition is the coalescence and 
percolation of randomly nucleated microcracks. 

In the analysis of the numerical results obtained by 
the computer simulation of our dynamical model, we re- 
ferred to the theory of percolation on finite lattices just as 
an example of critical phenomena occurring in finite sys- 
tems, without assuming percolation of microcracks. In 
this sense our analysis is free of model assumptions. In 
the snapshots of the final scenarios of collision processes 
it can be seen that the complicated elastic field arising 
due to propagation and interference of elastic waves gives 
rise to correlated crack growth. Hence, in our case the 
assumption of percolation of cracks cannot be valid, and 
the values of the exponents obtained numerically are dif- 
ferent from the corresponding exponents of percolation. 

An important feature of our model which can affect the 
result of a fragmentation process is the existence of ele- 
mentary, non- breakable polygons, which hinders us from 
observing fragmentation on lower scales. In Ref. us- 



ing a discrete model of solids similar to ours, detailed 
tests of the effect of the polygon size on the shape of the 
mass distribution of fragments in impact fragmentation 
were performed. It was found that if the tessellation of a 
solid is fine enough to reproduce the macroscopic elastic 
behavior correctly, the effect of the existence of unbreak- 
able polygons on the size distribution of fragments occurs 
solely in the vicinity of the cutoff size, and the overall 
shape of the distribution is not affected. This also im- 
plies that the phase transition behavior revealed by our 
simulations is not influenced globally by the element size, 
but it can be expected that the critical regime gets nar- 
rower when the element size decreases at a fixed value of 
the macroscopic size of the colliding bodies. 

Studying kinetic models of fragmentation phenomena 
in Refs. [|23|]3q | it was pointed out that under certain cir- 
cumstances a so-called shattering transition might occur, 
when a macroscopic amount of mass of the fragmenting 
system is transformed into a dusty phase, i.e. into single 
polygons in our case. Since in discrete dynamical models 
such a shattering transition can occur only at very high 
impact energies its investigation is not included in the 
present study. 

The binary break-up of bond percolation clusters due 
to the removal of a single bond has also been studied 
in the context of fragmentation phenomena. The char- 
acteristic quantities describing the fragmentation of per- 
colation clusters are the number of fragmenting bonds 
a s (p), i.e. the number of those bonds in the cluster of 
size s whose removal results in two disconnected clus- 
ters, and the probability P s / tS (p) that the break-up of a 
cluster of size s results in a daughter cluster of size s' 
p9[ . a s (p) and P s ^ s (p) are of fundamental interest be- 
cause they describe the connectivity of percolation clus- 
ters and they may also serve as an input for the rate equa- 
tions describing fragmentation processes as a sequence of 
binary events |2^,^|,^9|. At the critical point p = p c 
both a s (p c ) and P s ' jS (p c ) show interesting scaling behav- 
ior whose scaling exponents could be related to the criti- 
cal exponents of percolation |m| , but no phase transition 
was claimed in the break-up process. Although the form 
of P s '. s (p c ) as a function of the daughter mass s' and its 
dependence on the global size of the fragmenting object 
s [ |39| is similar up to some extent to the behavior of our 
fragment mass distribution at the critical point F(m;ri c ) 
(see Fig. ||), there is no direct analogy between the bi- 
nary break-up of percolation clusters and our dynamic 
fragmentation process. 

Up to our knowledge no systematic experimental study 
of the transition from damage to fragmentation in im- 
pact of solids has been performed so that we cannot con- 
front the results of the simulations with experiments. For 
a deeper understanding of the transition predicted by 
our simulations further experimental and more analyti- 
cal theoretical studies are needed. 
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TABLE I. The parameter values used in the simulations. 



Parameter 


Symbol 


Unit 


Value 


Failure elongation of a beam 


t e 


% 




3 


Failure bending of a beam 


t@ 


degree 




4 


Surface energy 


E s 


erg 


1.4 


• 10 6 


Damage threshold 


E d 


erg 


1.7 


■ 10 8 




Vo 


1 





.08 


Range of velocities 


Vo 


c m 

s 


200 


- 3500 



7 



TABLE II. Test of the consistency of the critical exponents. Starting from the three possible pairs of r, j3, and 7 the value 
of the other exponents were derived using the scaling laws. 

Cjyg) (t, 7 ) 08,7) 

0.11 ± 0.02 7-(r-2) ~ Q 09 ± q Q25 0.11 ±0.02 

7 (3 ~_2 ~ = °- 29 ± °- 08 0.26 ±0.02 0.26 ±0.02 

<r 2=2 =2.45 ±0.6 ^ =2.8 ±0.28 -^=2.7 ±0.21 

P 7 7+P 

r 2.27 ±0.05 2.27 ± 0.05 2 ± ^ = 2.3 ± 0.04 

a 2- fc^M = 1.48 ±0.12 2 - (t 3 ~!^ 7 = 1-54 ± 0.06 2 - (2/3 + 7) = 1.52 ± 0.045 
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FIG. 2. The final breaking scenarios of collisions of discs 
at different impact energies Eo- The radius of the discs is 
R — 20cm, while the average size of a randomly shaped poly- 
gon is 1cm. The lines connecting the center of mass of the 
neighboring polygons symbolize beams. The missing beams 
have already been broken. The values of the parameter r\ are 
0.09, 0.2, 0.3, 0.5 for a, b, c and d, respectively. 
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FIG. 5. The mass of the largest fragments divided by the 
total mass of the colliding system. The arrow indicates the 
fragmentation threshold rj c . 



7] 

FIG. 3. The energy released by beam breaking Er and the 
kinetic energy stored in the motion of fragments _Eo — Er. For 
comparison, the surface energy E s of the model, the binding 
energy of the sample Et and the energy of the impact Eo are 
also indicated. The energy unit is erg. For the exact values 
see Table @. 




V 

FIG. 4. The ratio m of the energy released by breaking Er 
and the total kinetic energy E . The transition point (frag- 
mentation threshold) between the damaged and fragmented 
states is identified with the position of the maximum of £r. 




T) 

FIG. 6. The energy release curve cr at different system 
sizes R. 
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m [g] 

FIG. 7. Mass distribution of fragments at different values 
of r] below and above the fragmentation threshold rj c . The 
peak of small fragments is due to the fact that single poly- 
gons are unbreakable in the model. The mass of the colliding 
discs with radius R — 20cm was 1250 gram. 
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FIG. 8. Mass distribution of fragments at the transition 
point rj c for three different values of the system size R. 
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FIG. 9. (Ml/M{} as a function of rj for five different values 
oiR. 
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FIG. 10. Scatter plot ol moments for the system of size 
R=20cm. The straight lines were fitted to the data points in 
the range M J 2 /M{ > 300. 




I Wc 1 

FIG. 11. Determination of the order parameter exponent 
13. The data of Fig. |^ is plotted here as a function of \q — r) c \ 
for 7] < T)c- The exponent /3 can be obtained as the slope of 
the least square fitted straight line. 
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FIG. 12. To obtain the exponent 7 we plotted < M 3 2 /M{ > 
belonging to the system size R — 20cm of Fig. ^ as a function 
of |?7 — ?7 C for rj > r/c. Data points are taken only from the 
vicinity of r\ c . The exponent 7 is obtained as the slope of the 
least square fitted straight line. 
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